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§1. Introduction 



We begin with briefly reviewing the conventional time-energy uncertainty relation in 
quantum mechanics. Let A be an ovservable without explicit time-dependence and \ip{t)) 
be a normalized quantum state vector obeying the Schrodinger equation with a hermitian 
Hamiltonian H. If we define AA and ta by 

AA = ^{m\A'\m) - mmm)' , a-i) 



Ta 

and take the equation 



, (1.2) 



into account, we are led to the uncertainty relation^) 

taAH > I . (1.4) 

The quantity ta is interpreted as the time necessary for the distribution of {tlj{t)\A\il;{t)) 
to be recognized to have clearly changed its shape. 

In contrast with the well known result given above, Anandan and Aharonov have 
recently succeeded in obtaining quite an interesting inequality.^) They consider the case 
that the I'lpit)) develops in time obeying 

{m\m) = ^, (1-6) 

where H{t) is an operator which is hermitian and might be time-dependent. They con- 
clude that 

/*' AS{t)dt > hATCcos{\{ilj{ti)\ijj{t2))\) , (1.7) 
Jti 

where AS{t) is given by 



As{t) = ^m)\Hmm)-m)\mm)? ■ (i-s) 

The inequality (1.7), which we refer to as the Anandan- Aharonov time-energy uncertainty 
relation, has been derived through a geometrical investigation of the set of normalized 
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quantum state vectors. The r.h.s. of (1.7) can be regarded as the distance between two 
points in a complex projective space. We note that Montgomery also proposed an inter- 
esting time-energy uncertainty relation of a similar nature. 

In this paper, we seek the generalized version of (1.7). We consider a set of N or- 
thonormal vectors Hipiit)) : i = 1,2, . . . , N} satisfying 

ii^iitMjit)) ^ 5ij , 1,3^ 1,2,..., N, (1.9) 

each of which obeying the Schrodinger equation (1.5). We define N xN matrices A{ti, 
and K{t-i, ^2) by 

A{ti,t2) ^ {aij{ti,t2)) , aij{ti,t2) ^ {A{tiM3{t2)) , (1.10) 

K{h,t2)^A^{tl,t2)A{h,t2) (1.11) 

and Kiiti, 12), i = 1,2, . . . ,N,io be the eigenvalues of K{ti,t2). Defining the generalization 
of (1.8) by 



we find that A£iv(^) satisfies 



N N 
\ i=l i,j=l 



[ ' AEN{t)dt >h.Yl |Arccosyt(i]7^2)| . (1.13) 
The inequality (1.13) can be written in an operator form as 

r ^Tr{P{t)[H{t),[H{t),P{tmdt 

(1.14) 

> V2h^Tr{{Arccos^P{h)P{t2^) , 

where P{t) is defined by 

N 

p{t) = j2m)){mi (1-15) 

i=l 

and Tr denotes the trace in the Hilbert space. The result (1.13) is obtained through a 
geometrical investigation of the Grassmann manifold mentioned below. 

This paper is organized as follows. In §2, we introduce some objects such as geodesic, 
distance, etc., defined on the set, Gn, of A'"-dimensional linear subspaces of a Hilbert 



space. In §3, we obtain an explicit formula for the distance between two points in G^v- 
In §4, we discuss that the distance introduced above satisfies the properties of distance 
including the triangle inequality. The inequality (1.13) is derived in §5. The final section, 
§6, is devoted to discussions. Some appendices are attached to explain some necessitated 
calculations. 
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§2. General discussions of the distance in Gn 
2.1 A'^-th Grassmannian 



Given a Hilbert space h,we consider vectors i = 1,2,..., A/", belonging to h and 
satisfying {ipi\ipj) = Sij. We call the set 

^r = (|V'i),|V2),...,|V'iv)) (2.1) 

an A^-frame of h and the set 

[qj] ^{-qju-.ue U{N)} (2.2) 
an iV-plane of h,where "^u is defined by 

TV AT N 

*M = (XI J2 I^j>i2, \^k)ukN)- (2.3) 

1=1 j=l k=l 

It is clear that the [\E'] and the projection operator 

TV 

p=J:\i^^)m (2.4) 

i=l 

are invariant under the replacement ^ — > ^u. We denote the set of all the ^'s of h by 
Sn- Then the set Gn defined by 

Gn = m -.^eSN} (2.5) 

is known to constitute a manifold of complex dimension A^(dim h — A^). We hereafter call 
Gjv the A"-th Grassmann manifold, or simply the A"-th Grassmannian. 

2.2 Geodesies in the set of unitciry operators on h 

We denote the set of unitary operators on hhj Q. A local coordinate of Q is denoted 
by s = (s^, s^, s^, . . .), s^, s^, s^, . . . G R. We define the infinitesimal distance between two 
unitary operators W and W + dW by 

D{W,W + dW) = \\dW\\, (2.6) 
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where the Hilbert-Schmidt norm \\w\\ of an operator w is defined by 



\w\ 



Tr denoting the trace on h. If we define g^uis) by 



d 



(2.7) 



(2.^ 



we have 



D(W, W + dW) = ,Jg^J^s)ds^^ . (2.9) 

The g^iyis) defined above can be regarded as the metric tensor of the set Q of unitary 

operators on h. It is evident that g^u{s) transforms as a tensor under the transformation of 
the coordinate s and the value of the infinitesimal dcstancc D{W, W + dW) is independent 
of the choice of s. Given the metric g^y{s), a geodesic of VL is defined as a solution 
s{t) — {s^{t),s^{t), . . .),t eR, of the equation 



s^it) + (5(t))s^(t)s<^(t) =0, = 1, 2, . . . , 
rL(«) = l9'''{s){dxg..{s) + d,gx.{s) - d.gx,{s)}, 



(2.10) 
(2.11) 



where {g'^^{s)) is the inverse of {g^„{s)) and the dot denotes the derivative with respect 
to t. 

The one parameter subgroup {e**^ :Y — Y\tE R} should be more or less related to 
a geodesic of Q since we have 



dW{t) 



dt 



and 



d 



itY 



dt 



l{iim=o. 



(2.12) 



(2.13) 



The fact is that any geodesic in Vt passing the point 1 (unit operator) can be regarded 
as a one parameter subgroup of the above form. Although this fact can be seen in 
mathematical liter atures,^) we discuss it in the Appendix A for self-containedness. We 
note that the length of the geodesic {e**^ : Q <t <1,Y — Y^} connecting the two points 
1 and e'^ in Q, is given by 

" dt^ \\Y\\. (2.14) 







Jo 


dt 
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2.3 Distance in Gat 

To an iV-frame ^(t) = \ip2{t)), • • • , |^7v(t))) e S'^v, < t < 1, there correspond 

an A'"-plane [^{t)] G Gn and a projection operator P(i) = J2iLi Since the 

eigenvalues of P(l) are equal to those of P(0) including multiphcities, there exist a unitary 
operator W such that 

P{1) ^W^P{0)W, W^e^^, ^Y. (2.15) 

The discussion of the previous subsection suggests that we might be able to define the 
distance d([^'(0)], [^1/(1)]) between two points [^1/(0)] and [^(1)] of the Grassmannian Gn 
by the formula similar to the r.h.s. of (2.14). We define it by 

d([*(0)],[*(l)]) = Min||F|| (2.16) 

where E is the set of hermitian operators specified by -P(O) and -P(l) in the following way: 

E = {y : y = Y{P{0), P(l)) = -Y{P{1), P(0)) = e-'^P{0)e^^ = P(l)}. (2.17) 

Supposing that Min ||y||, Y e E, is attained hy Yq e T,, we have 

^([^(0)], [*(!)]) = roll- (2.18) 

As will be shown later, if we require that the functional form of yo(-P(0), P(l)) should 
be fixed independently of the choice of P(0) and P(l), the yo(-P(0), P(l)) is determined 
uniquely. It should be stressed that the r.h.s. of (2.16) or (2.18) is invariant under the 
replacement — > G U{N) and can be regarded as a quantity defined on Gn- After 

obtaining the explicit expressions of Yq and o?([^(0)], [^(1)]) in §3, we will discuss in §4 
that the above defined distance in Gn satisfies the property of distance: 

dm,m) = dm,m>o, (2.19) 

d([*],[$]) = o ^ m = m, (2.20) 
dm,m)<dm,[E])+dm,m, (2.21) 

for any G Gn- 
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§3. Explicit formula for the distance in Gjv 

Adopting abbreviated notations Pq and Pi for the projection operators P(0) and -P(l), 
respectively, we have 

Pl^Pm, (Pmy^Pm, m = 0, 1. (3.1) 

It is not difficult to see that the most general form of Y satisfying Y = Y{Pq, Pi) = 
-Y{Pi, Po) = Y^ is given by 

Y = a{PoPi)Po - a{PiPo)Pi + i{l3{PoPi)P,Pi - /5(PiPo) A^o}, (3.2) 

where a{z) and (5{z) are real analytic functions of z involving no inverse powers of z. As 
is explained in the Appendix B, we obtain 

; X 2 

N 



IV(e^^Pie-^Po) ^Y.^i{ cos - (1 - ^i)/3(^i) / , (3.3) 

where Ki,i — 1,2, ... ,N, are the eigenvalues of the N x N matrix K whose ij-element is 
given by 

TV 

fe=i 

The s{k,) in (3.3) is calculated from (B.23) and (B.25) to be 

s{k) = (1 - K){{a{n)y + K{P{K)y). (3.5) 
In the Appendix C, the eigenvalue Ki is shown to satisfy 

0<Ki<l, i^l,2,...,N. (3.6) 
Another requirement e~*^Poe*^ = Pi characterizing F G S is equivalent to 

Tr(e^^Pie-^^Po) = N (3.7) 

as is seen from 

\\e-'^Poe^^ - Pif = 2{N - Tr(e^^Pie-*^Po)}. (3.8) 

Prom (3.3) and (3.7), we obtain the condition to specify the functions a{z) and (3{z) 
defining y e E: 

N 



Y^C{Ki) cos\,Js{Ki) - <P{Ki)) = N, (3.9) 



1=1 



where C{k) and ^(k) are given by 

^^^^ - " {a{K)Y + k{(5{k)Y ' ^^-^^^ 
i^nm = . (3.11) 

Except for the trivial case = K2 = ■ ■ ■ = kn = 1, the condition (3.9) is reahzed only 
when a{Ki) = 0,i = 1,2, . . . , N. Since we are considering generic cases where /tj's satisfy 
(3.6) and requiring that the functional form of a{z) is independent of Pq and Pi, we 
conclude that a{z) vanishes identically : 

a{z) = 0. (3.12) 

Then C{k) equals 1 and we see from (3.5), (3.9) and (3.11) that P{k) should satisfy 



tsin(p{K) = -J- — -sgn{(3{K)) (3.13) 

V K 

and 

-{</.(«:) - ^JK{l-KmK)\} e Z. (3.14) 
Eliminating (j){K) from (3.13) and (3.14), we are led to the condition 



taii(v/41-^)/3W) = -/^ 



' : ' (3.15) 

to determine P{z). We have seen that the operator Y belonging to E of (2.17) is specified 
by (3.2), (3.12) and (3.15). 

We next determine Yq which attains the minimum of e E. It is now easy to 

observe that Yq is uniquely given by 

Yo = i{Po{PoPi)PoPi - Po{PiPo)PiPo}, 
= GiPo,P,) 

where Po{z) is given by 



(3.16) 



Poiz) = , Arctan\/ 



A.iccos\fz (3-17) 



^_-_(2n)!^ 



With the help the relation (Pq - Pif^'PoPi = (1 - PqPiYPqPi, n = 0, 1, 2, . . . , we have 
G{P., P^) = ]f: J^fL ^Po - PiriPoPi - PiPo). (3.18) 

The distance (i([^(0)], [^(1)]) is calculated by (2.16). Through procedures similar to those 
of Appendix B, we obtain 



d{[^m, [^(1)]) = V2Tr((a: - x^){Po{x)y), x = PoPi, 



N (3.19) 

. 2^(Arccosv^)^. 
\ i=i 

In the case of = 1, the above distance reduces to \/2Arccos|(-?/'(0)|-?/'(l))| and repro- 
duces the distance of the complex projective space, which was utilized by Anandan and 
Aharonov.^^ After some manipulations, we obtain an important relation 

||[G(Po,Pi),Po]|| = ||G(Po,Pi)||, (3.20) 

which will be useful for later discussions. 

In a recent paper, Avron, Seller and Simon^^ argued algebraic properties of a pair of 
projection operators. They noted that, under a suitable condition on the norm of Pq — Pi, 
the unitary operator 

^1 - (P„ - F,)' 

transforms Pq iiito Pi '■ 

U^PqU = Pi. (3.22) 

The expression on the r.h.s. of (3.21) makes sense since (Pq — PiY commutes with Pq and 
Pi. It can be read off from Ref. 6) that (3.21) was originally obtained by T. Kato many 
years ago. In the Appendix D, it is shown that our exp(iG'(Po, Pi)) is nothing but the U 
of (3.21) : 

giG(Po,Pi) _ u_ ^3 23) 
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§.4 Property of distance 



We here discuss how the properties (2.19), (2.20) and (2.21) are assured. It is evident 
that the d([*], [$]) defined by (2.16) staisfies (2.19). Noticing that, in the case that 
[if] — [$], all the eigenvalues corresponponding to KiS of §3 are equal to 1 and all the 
Arccosy^, i = 1, 2, . . . , A/", vanish, we see that (2.20) is also satisfied. 

We now go on to the triangle inequality (2.21). For the case of = 1, (2.21) is 
equivalent to 

Arccosy^ < Arccos\/a + Arccos\/A, (4.1) 
where k, X and a are defined by 

K^mi^)f, a^mof, A=i(ei0)r (4.2) 

with 10), lip) and |^) being unit vectors. The inequality (4.1) can be proved algebraically 
as is seen in Appendix E. In general cases of A^ > 2, we have not succeeded in proving 
(2.21) algebraically relying solely upon the formula (3.19). In the following, we describe 
an analytic proof of (2.21). 

We consider three A^-planes [^], [S], [$] G Gat and denote the projection operators 
corresponding to them by P*,Ps and P$, respectively (P^ = Sili etc.). We 

define hermitian operators Yi, Y2 and by 

Y, = -G{P^, P^), Y2 = -G{P^, Ps), n = -G{Ps, P$), (4.3) 

where G{Po, Pi) is defined by (3.16). A piecewise smooth geodesic T = {^{t) : < t < 1} 
in fl is defined as 

exp (i-Y2) ; 0<t<to 

^^^) = ^ /t-i X ^4.4) 
exp (i- — -Ys ] e^^'; to<t<l, 



I -to J 

where to is a constant satisfying < to < 1- 
We define a set of projection operators by 

Ajv = {P : = P, P^ = P, TrP = A^} (4.5) 

and consider a path 

C = {P{t) : < t < 1} C Ajv, (4.6) 
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where P{t) is given by 



Pit) = ^it)p^^ity. 



(4.7) 



The projection operator P{t) satisfies P{0) = Pv&, P(to) = Ps and P{1) = P^- The length 
of C is defined by 



1{C) ^ [ \\dP\\ ^ [ 
Jc Jo 



f 


dP{t) 
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dt 



dt. 



(4.8) 



Then we have 



m = rU[Y,,p^]\\dt+ f ^\\[Y,,p^]\\dt 

Jo to Jto i- — to 

^\m,p^]\\ + \m,p^]\\ 
= 11^211 + iimi 

= dm,m+dm,m), 

where the relation (3.20) has been made use of. Similarly, the length of the path 



(4.9) 



Co = {Po{t) : < t < 1} C Aiv, 



with 

is given by 



Co 



\\dP\\ 



Po{t) 



dPoit) 



JtYi 



Pis,e 



-itYi 



(4.10) 



(4.11) 



dt 



dt = \\[Y,,p^]\\ = \\Y4 = dm,m). (4.12) 

The metric tensor G^u{s) corresponding to an infinitesimal distance \\dP\\ is given by 

\\dP\\ = ^/G^Js)ds^, (4.13) 



G^,{s)^Tr{d^P{s)d,P{s)), 



d 



(4.14) 



where s denotes a set of local coordinates of A^. The /(Co) is stationary under any 
variation 6 of Cq with the end points P^ and P$ fixed. This is because the length /(Co) is 
a function only of Pq, and P$. On the other hand, for some choices of P^{e.g., — 
{ii\(t>j) = 0, i, j = 1, 2, . . . , N), the /(C) can be clearly larger than /(Cq). For generic P=, 
the /(C) is not stationary under the above mentioned variations 5 since the value of /(C) 
varies continuously with P^. For some choices of Ps(e.g'., Ps = P$ or P3 = P^), the /(C) 
coincides with /(Co). From these discussions, we see /(C) > /(Cq), implying (2.21). 
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Note that the Q{t) defined by 



Q{t) = tP^ + (1 - t)P9, 0<t<l, 



satisfies 



dt 



dQ{t) 
dt 



dQit) 
dt 



dt 



(4.15) 



= and 6{1{Q)) = with 1{Q) given by 1{Q) = /q 
P^W < /(Co). The set {Q{t) : < t < 1}, however, does not belong to Ajv since {Q{t)y 
is different from Q{t) for general values of t. So, we should not regard ||P$ — P^\\ as the 
distance in A^r between P^ and P$ although we have 1 1 P$ — 1 1 < 1 1 P$ — Ps 1 1 + 1 1 Ps — P* 1 1 ■ 
Finally, we mention one more inequality. If we define a hermitian operator Y{t) by 



7(i)=e'W*), 0<t<l, 

R{t) eU, R{0) = 0, i?(l) = 1, 
ITOII = const. = (1)11, 
the relation e'^^e'^^ = 7(1) = e*^*^^'' and the Campbell-Hausdorff formula yield 

tY{l) = tYs + 1Y2 + ^[F3, Y^] + ^([[Fs, F2] - [[V3, ^2], n]) + • • • • 



(4.16) 

(4.17) 
(4.18) 

(4.19) 



The path Ci = {Pi{t) : Pi{t) = e'^^^^^P^e''^^^^^ also connects P^ with P$ in Ajv- 
The Pi(t) satisfies | 



dPijt) 
dt 



= which is necessary but not sufficient for Ci to be a 
geodesic. The /(Ci) for generic Ps, however, would not be stationary under the variations 
6 mentioned above. Discussions similar to those of the case of 1{C) lead us to ||y(l)|| = 
/(Ci) > /(Co) = 11^1 II- On the other hand, as is seen in the Appendix F, we have 



dj{t) 



dt 



> 



dR{t) 



dt 



\Y(1] 



(4.20) 



From (4.8), (4.9), 



dP(t) 




dl{t) 


dt 




dt 



\Y,, 



implied by (3.20), (4.17) and ||r(l)|| > ll^iH, we obtain 

+ iinii > r(i)ii > rill, (4.21) 



which somewhat reveals more detailed feature of ||y2|| + H^all > H^il 
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§5. Time-energy uncertainty relation 



Having obtained the explicit formula (3.19) for the distance in Gjv, we here derive the 
uncertainty relation (1.13) or (1.14). The projection operator P{t) is defined by (1.15) 
and \tlji{t)),i — 1,2, . . . ,N, develops in time obeying (1.5). We then have 

P(t-rdt) = P{t)-r'^[H{t),P{t)] 

When [^(0)] and [^(1)] in (3.19) are close to each other, Ki,i — 1,2,...,N, are nearly 
equal to 1. Noticing that (Arccos-0i)^ 1 — k for k f=i 1, we see 



dim)], m + m 



N 

2^{1 -K,{t)), (5.2) 



i=l 



where Kj(i)'s are obtained from P{t) and P{t + dt) by similar procedures to those of 
previous sections. Since, in the above case, we have TrP{t) — N and 

N 

Tr{P{t)P{t + dt))^J2'^iit)^ (5-3) 
1=1 

(5.2) can be rewritten as 

d{m)],mt + dt)]) = ^j2TriP{t){Pit)-Pit + dt)}). (5.4) 
Eqs. (5.1), (5.4) and the relation Tr([A, B]) = yield 

dim)], mt + dt)]) = ^-f^Tr{P{t)[H{t),[H{t),P{t)]]) 

^\^TT{[P{t),HmH{t),P{t)]) ^^^^ 



dt 

= \\dpm- 

It can be easily seen that the r.h.s. of (5.5) is proportional to AS^^t) defined by (1.12). 
Now we are led to 

d{[^{t)],[^{t + dt)]) = ^ASN{t)\dt\. (5.6) 
For finitely separated [^(^i)] and [^(^2)] in Gn, the triangle inequality (2.21) imphes 

j^*' ASj,{t)dt > -^d([*(ti)], [^{t2)]),t2 > ti. (5.7) 
The formula (3.19) then leads us to (1.13) or (1.14) or (1.15). 
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§6. Discussions 

In this paper, we have been mainly engaged in obtaining the formula (3.19) for the 
distance in Gjv. As an application of it, we obtained the generalized version of the time- 
energy uncertainty relation, (1.13), of Anandan-Aharonov type. For A'"=l, it reduces to 
the result of Anandan and Aharonov, (1.7). 

Our definition of the distance c?([^I/], [$]) in Gjvj (2.16), is intimately related to that 
in Q. The discussions around (4.13) or (5.5) indicate that the distance between two 
infinitesimally separated A^-planes [if] and + d'if] is given by 

d([*],[^ + d^]) = lldP^II, (6.1) 

dP-^ — — P-qt, (6-2) 

where P-^+d* and P<i( are projection operators associated with the //-frames \E' + d"^ and 
"if, respectively. If we make use of some real coordinates s — (s^, s^, . . .) to specify [^] or 
P^, we have 

{dm, + d*])}' = G^,{s)ds^ds\ (6.3) 

with G^v{s) given by (4.14). Metric tensors of this kind have been discussed recently by 
two of the present authors.^) We stress that the metric tensor and the geodesic for the 
Grassmannian can be simply expressed in terms of projection operators belonging to 

^N. 

When the Hamiltonian H is independent of time, the Afjv(t) in (1.12) does not depend 
on t and is determined only through [^(0)]. We suppose that '^{t) is orthogonal to 
*(0), ie., (V'i(i)lV'j(O)) = 0, i, j = 1, 2, . . . , and that A£:iv(0) is nonvanishing. Then we 
have = 0, i = 1, 2, . . . , A", and hence 

iA£:Ar(0) > 7rV]V^. (6.5) 

In other words, an A"-frame ^ needs 'K\/Nh/ {2AEn{^)) of time to develop to the one 
orthogonal to the original one. 

Other applications of the distance formula (3.19) will be discussed elsewhere. 
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Appendix A. Geodesic and one peiranieter subgroup in Q 



(a) In the following,we discuss that a geodesic in Q, 

r = {7(i) : t e 11,7(0) = l(unit operator)}, 



(^.1) 



can be identified with a subgroup {e**^ : Y = Y\ t € H} of ^2. More systematic expositions 
can be found in Ref. 5. 

(b) For any uj £ fl,we define the mapping as follows. For arbitrary a;, a;' e fl, we can 
think of a geodesic F such that u; — 7(a) and to' — 7(6), a, 6 e R. We define i^^ by 



lMt))=^{2a-t). 



{A.2) 



Eqs. (2.10) and (2.12) imply that ||^|| is constant on F. Then ,we can characterize the 
mapping /^^ by 



and 



|U7(^))IH 



dt 



,i eR. 



It is clear that Iuj{uj') is well-defined for any uj' E Q and that is unique, 
(c) On the other hand,if we define the mapping J^^ by 



we have 



d 
It 



{ d 



t=a 



dt 



uj{'y{t))-'Lj) 



t=a 



-M(t)-^^7(^)-^a;},.„ 



dt 



,dl{t) 
' dt 



)t=a 



and 



Ui{t)) 



-^^{tr^^itru; 





dl{t) 




dt 



,i eR. 



(^.3) 
{AA) 

(A.5) 



(A.6) 



(^.7) 



{A.i 



{A.9) 
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Thus, we can identify J^; with 

Jo, = /a,- (AlO) 

(d) Now, we assume that the unit operator 1 hes on F and choose the parameter t such 
that 7(0) = 1. We have 

A(*)(A(o)(7M)) = I,itM-u)) = ^{2t + u) (All) 

and 

^7{t)(-^7(o)(7(^^))) = J,(tMu)-') = ^(thiuHt). {A.12) 
From (A. 10-12), we obtain 

7(2t + «)=7(t)7(«)7(t) (A13) 
Putting u — 0,t, 2t, 3i, . . . , in (A. 13), we have 

7(nt) = {7(t)r,n = 0,l,2,.... (A14) 
Since 7(i) should be continuous in t, we conclude that j{t) can be written as 

7(t) = e**^, ten, Y = YK {A.15) 
We thus arrive at the statement in (a). 



Appendix B. Calculation of Tr(e*^P(l)e-*^P(0)) 

In this Appendix,we calculate Tr(e*-^P(l)e~*^P(0)) through several steps, where Y is 
the operator given by (3.2). We use the abbreviated notations Pq = P(0), Pi = P(l), a; = 
P(0)P(1),|/ = P(1)P(0). 
(a) Making use of the relations 

xPi — X, PqX — X, 

(P.l) 

Z/"Pi = Prx^, y^x = Pix", n = 1, 2, 3, . . . , 
repeatedly, we find that Y in (3.2) satisfies 

Y'^Pi = fn{x) + Pignix),n = 0, 1, 2, ... , {B.2) 
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Hence, we have 



with fn(x) and gn(x) given by 

a{x) + i(3{x)x {a{x) + if3{x)}x \ 
-{a{x) + il3{x)) —{a{x)+il3{x)x) J' 

e^^Pi = (1 Pi)e^«(^) 5 j . 
Similarly, we obtain 

= f„{y) + P„g„{y), 

(|j^j)M«wr(?).»^o,i,2,..., 

-{a{y) + tP{y)y) -{a{y) + tP{y))y \ 

(^{y) + iP{y) (^{y) + iP{y)y /' 



Q{y) 

and 



e-^^Po = (1 Po)e-^'3fe) ( J 



(b) From (B.2) and (B.6),we obtain 



B.3) 
BA) 

P.5) 

P.6) 
P.7) 

P.8) 
B.9) 



Tr(e'^Pie-^^Po) = ^ ^—^c^^, (i^.lO) 
where Gmn is given by 

= 'IV({/„(x) + Pi^„.(x)}{/„(y) + PoUy)})- (B.n) 
Noticing that (B.3) and (B.7) ensure the existence of hm{x) and hn{y) such that 

fm{x) = xhm{x), fn{y) = yK{y),n = 0, 1, 2, ... , (-B.12) 

we obtain 

Gmn = +5n(y)})- (5-13) 

In deriving (B. 13), the following relations have been made use of: 

Tr{xP+^y'^+^) = TT{PixPy''+^) = Trix^+^Poy") 

= TriPixPPoy'^) = Tr{yP+i+^) = Tr(xf+«+^), (P.14) 
= 0,1,2, 

18 



(c) Prom the definitions of Po,Pi and y,we have 
We then have 

AT 

y'^'- E \A{l)){AK%{^j{0)\,q = 0,1,2,..., {B.17) 
where A is the N x N matrix whose ij-element is Uij and the N x N matrix K is defined 

by 

K = A^A. {B.18) 

Prom (B.17), we obtain 

IV(y^)=tr(X«),g = 1,2,3,..., (R19) 

where tr denotes the trace of an iV x iV numerical matrix. Eqs. (B.10),(B.13) and (B.19) 
yield 

Tr(e*^Pie-^^Po) = tr{KG{K)G{K)), (5.20) 
where G{K) and G{K) are given by 

GiK) = E ^(^n(i^) + 9n{K)) = (X-^ l)e^«(^) ^ , (S.21a) 

n=0 V / 

G{K) = £ ^^(/in(i^) + UK)) = {K-' l)e-^(^) ( ° V {B.2lh) 

n=0 V / 

The apperrance of K^^ in (B.21) is only spurious and no dangerous procedure is involved 
there. 

(d) If we denote the A^-th unit matrix by 1n,Q{K) can be written as 

Q{K) ^Y.q,{K)a,, {B.22) 

1=1 

q^{K) = ]^{K-U){a{K) + i(5{K)), 

q,{K) = '-{K + U){a{K) + iP{K)), (^-23) 
q3{K)^a{K)+iKP{K). 
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We now have 



G{K) = cos Js{K) + {K- 1n)P{K)- 



sin s{K) 



with s{K) given by 

s{K) = j:{q^{K)f 

Similar manipulations yield 

G{K) = G{K). 
From (B.24),(B.26) and (B.20),we obtain 

/ 



Tr(e^'^Pie-^^Po) = tr 
yielding (3.3). 



I f 

K I cos ^s(X) + {K- In)I3{K)- 



sin J s{K) 



\ 



s{K) 



(S.24) 



(S.25) 



(S.26) 



{B.27) 



Appendix C. Proof of (3.6) 

It is clear that K is hermitian and positive definite so that we have Hi > Q,i — 
1,2, . . . , N. Supposing that uKu^ ,u G U{N),is diagonal,we obtain 

N 

Uikihi - Kki){u^)ij = (1 - = (^kj), (CI) 

k,l=l 

where |q) is given by 

N 

Iq) = EK)*I^;), (C-2) 

1=1 

k=l 

Putting i = j in (C.l),we have Ki < l,i = 1,2, . . . , N. 
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Appendix D. e^G{Po,Pi) ^nd U 
Prom (3.18), we have 



OO -1 

,iG{Po,Pi) ^ ^ 



n=0 ^■ 



L = PoPi - PiPo = -(Po - Pi)(l - n - Pi) p.3) 
since J commutes with L. By making use of the relations^) 

{Po-Pi,l-Po-Pi} = 0, {DA) 
(Po - Pif + (1 - Po - Pi)' = 1 {D.h) 

and 

^1 - (P„ - Pi)2 

we are led to 

1 - Po - Pi 



giG(Po,Pi) ^ cos(Arcsin(Po - Pi)) - sin(Arcsin(Po - Pi)) 
It is easy to derive (3.21) from (D.7). 



^1 - (Po - Pi) 



2 
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Appendix E. Proof of (4.1) 

The definitions (4.2) of a and A mean that and |0) can be expressed as 

10) = VAe^^lO + \r]) {E.2) 

where 5, e e R and k) and I77) are appropriate vectors orthogonal to |^). Since and 
10) are unit vectors,we have 

(^k) = l-(7,(?7|7?) = l-A (£;.3) 

and 

|(r?k)| < vT^yr^. (£;.4) 

The lower bound of ^/k is estimated as 

= |(VAe-(e| + (77|)(v^e^^|0 + k))l 

= |^/A^e''^-^^ + (r/k)| (^-5) 

>v^-|(r/k)| 

> VAct - Vl - crVl - A 

Noticing the formula Arccos(-\/Aa- — Vl — A-\/l — cr) = Arccos-\/a + Arccos\/A, we obtain 
(4.1). 
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Appendix F. Proof of (4.20) 



We consider a two-parameter family of unitary operators 

W^(r,t) = e-^W = (s^(r,t,),s2(r,t, ),...), 

where s'^(r, t), = 1, 2, . . ., are local coordinates of Q and Y{t) is so chosen that 

d 



We note that (F.2) is equivalent to 



d_ 
di 



dW{r, t) 



(F.l) 



(F.2) 



dr 



(F.3) 



dst^ir, t) ds^ir, t) 



(FA) 



Given g^iy (s) by (2.8), we define A{r,t) by 

Air,t) ^ g,.is) ^ 
We have 

^(0,t)=0 (F.5) 

since W{Q,t) — 1 and s'^(0, i) is independent of t. On the other hand, with the help of 
the formula 9«g'^^(s) = Ti'^^{s)gpu{s) + T''^^{s)gpi,{s), we obtain 



-Air, t) = g,.{s) -^K^ir, t) + -^L^{r, t) 



where K^{r,t) and V{r,t) are given by 



K^(r f\ _ dh'^{r,t) ds^{r,t)ds\r,t) 



Qj.2 



dr dr 

ds'^{r,t) ds^{r,t) 



(F.6) 



(F.7) 



(F.8) 



drdt ' ^"^^"^ dr dt 
The fact that {VF(r, t) : r e R, fixed} is a geodesic in Q implies K''{r, t) — 0, while the 
assumption (F.2) or (F.3) means that 



d'tV'''^'^~d~r d^ ~ 3^^^')^^^ ^^'^^ - °- 



We now have 



d_ 
dr 



A{r,t) = 0. 



(F.9) 



(F.IO) 
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Eqs. (F.5) and (F.IO) lead us to 



A(r, t) = 0. 



(F.ll) 



For the 7(t) in (4.5), we obtain 



d-f{t) _ f dW{r,t) ^ dR{t) ^ f dW{r,t) 



dt 



dr 



r=R(t) 



dt \ dt 



r=R(t) 



From {F.ll) and (F.12), we see that 



(F.12) 



dl{t) 



dt 



dR{t) 



dt 



dW{r, t) 



dr 



+ 



r=R{t) 



dW{r, t) 



dt 



r=R(t) 



(F.13) 



We conclude (4.20) since we know 



dW{r,t) 
dr 



\Y{t)\\ = \\Y{1)\\. We note that a more 



general version of the above discussion can be found in Ref. 5. 
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